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Abstract 

We calculate the electric dipole form factor of the nucleon that arises as a low-energy 
manifestation of time-reversal violation in quark-gluon interactions of effective dimension 
6: the quark electric and chromoelectric dipole moments, and the gluon chromoelectric 
dipole moment. We use the framework of two- flavor chiral perturbation theory to one 
loop. 
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Electric dipole moments (EDMs) [TJ |2] provide stringent bounds on sources of 
time-reversal (T) violation beyond the phase of the quark-mixing matrix [3J. Ex- 
periments are in preparation [4] which aim to improve the current bound on the 
neutron EDM, \d n \ < 2.9 ■ 10~ 26 e cm [5], by nearly two orders of magnitude. Novel 
ideas exist [6] also for the measurement of EDMs of charged particles in storage 
rings, including the proton — for which an indirect bound, \d p \ < 7.9 ■ 10~ 25 e cm, 
has been extracted from the atomic Hg EDM [7] — and the deuteron. Since the 
Standard Model prediction [HI [9] is orders of magnitude away from current experi- 
mental limits, a signal in this new crop of experiments would be a clear sign of new 
physics. 

The momentum dependence of an EDM is the electric dipole form factor (EDFF). 
Together with the well-known parity (P) and T-preserving electric and magnetic 
form factors and the P-violating, T-preserving anapole form factor, the P- and T- 
violating EDFF completely specifies the Lorentz-invariant electromagnetic current 
of a particle with spin 1/2. Although the full momentum dependence of a nuclear 
EDFF will not be measured anytime soon, the radius of the form factor provides 
a contribution to the Schiff moment (SM) of the corresponding atom, because it 
produces a short-range electron-nucleus interaction. 

There has been some recent interest [TUl - fT2"] on the nucleon EDFF stemming from 
the lowest- dimension T violation in strong interactions, the QCD 9 term. As other 
low-energy observables, both the EDM and the SM of hadrons and nuclei are difficult 
to calculate directly in QCD. Attempts have been made to extract the nucleon EDM 
from lattice simulations [13], but a signal with dynamical quarks remains elusive. 
One possible way to extract the EDM in this case relies on a extrapolation of the 
EDFF to zero momentum, which provides another motivation to look at the EDFF. 
QCD- inspired models have also been brought to bear on the nucleon EDFF [TTj . 

We would like to use a framework flexible enough to formulate the nucleon EDFF 
in the wider context of other low-energy T-violating observables such as the EDMs 
of nuclei. Such framework exists in the form of an effective field theory, chiral 
perturbation theory (ChPT) [HHT6]. (For introductions, see for example Refs. [171 
US].) Since it correctly incorporates the approximate chiral symmetry of QCD, 
ChPT provides not only a model-independent description of low-energy physics but 
also the quark-mass dependence of observables, which is useful in the extrapolation 
of lattice results to realistic values of the pion mass. The nucleon EDFF from the 9 
term has in fact been calculated in this framework [TU| [T2"] . and some implications 
of the particular way the 9 term breaks chiral symmetry were discussed in Ref. [19J. 
(For earlier work on the neutron EDM in ChPT, see for example Refs. J2UJ l2"Tj.) 
The momentum dependence of the EDFF is given by the pion cloud pUl [22]: the 
scale for momentum variation is the pion mass and the SM is determined by a 
T-violating pion-nucleon coupling. Assuming naturalness of ChPT's low-energy 
constants (LECs), one can use an estimate of this coupling based on SU(3) symmetry 
to derive [20] a bound on 9, 9 < 2.5 ■ 1CT 10 [12] from the current limit on the neutron 
EDM. ChPT extrapolation formulas for the nucleon EDM in lattice QCD can be 
found in Ref. [23]. 

The smallness of 9 leaves room for other sources of T violation in the strong inter- 
actions. Here we calculate in ChPT the nucleon EDFF arising from the effectively 
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dimension-6 interactions involving quark and gluon fields that violate T 
the quark electric dipole moment (qEDM), which couples quarks and photons; the 
quark chromoelectric dipole moment (qCEDM), which couples quarks and gluons; 
and the Weinberg operator, which couples three gluons and can be identified as the 
gluon chromoelectric dipole moment (gCEDM). These higher-dimension operators 
can have their origin in an ultraviolet-complete theory at a high-energy scale, such 
as, for example, supersymmetric extensions of the Standard Model. We construct 
the interactions among nucleons, pions and photons that stem from the underlying 
quark-gluon operators and use them to calculate the EDFF to the order where the 
momentum dependence first appears. As we will see, the sizes of the proton and neu- 
tron EDMs and SMs partially reflect the underlying sources of T violation. While 
much effort has already been put into estimating the EDMs from these sources [U [2] , 
the full EDFF apparently has been previously considered only within a particular 
chiral quark model [26]. Other implications of the different chiral transformation 
properties [27] of the dimension-6 operators will be considered in a subsequent paper 

EE]. 

Well below the scale My characteristic of T violation, we expect T- violating effects 
to be captured by the lowest-dimension interactions among Standard Model fields 
that respect the theory's SU(3) C x SU{2) L x U(l)y gauge symmetry. Just above 
the characteristic QCD scale Mqcd ~ 1 GeV, strong interactions are described by 
the most general Lagrangian with Lorentz, and color and electromagnetic gauge 
invariance among the lightest quarks (q = {u d) T ), gluons (G^), and photons (A^). 
The effectively dimension-6 T-violating terms at this scale can be written as 

C f = --q (d Q + d 3 r 3 ) a^ 5 q F, y - % -q (d + d 3 r 3 ) a^ 5 \ a q G% 



i dw tiv\a rabc^ia sib,psic 



in terms of the photon and gluon field strengths and G a , respectively, the 
standard products of gamma matrices 7 5 and in spin space, the totally antisym- 
metric symbol e^ vXa , the Pauli matrix r 3 in isospin space, the Gell-Mann matrices 
A a in color space, and the associated Gell-Mann coefficients f abc . In Eq. ([!]) the 
first (second) term represents the isoscalar d (d ) and isovector d 3 (d 3 ) components 
of the qEDM (qCEDM). Although these interactions have canonical dimension 5, 
they originate just above the Standard Model scale My/ from dimension-6 operators 
[23] involving in addition the carrier of electroweak symmetry breaking (the Higgs 
field). They are thus proportional to the vacuum expectation value of the Higgs 
field, which we can trade for the ratio of the quark mass to Yukawa coupling, m q / f q . 
Writing the proportionality constant as e5 q f q /Mp (47r5 q f q /Mj,), 



di~Oi e6—, , di ~ 4tt5— . , (2) 





in terms of the average light-quark mass m and the dimensionless factors 5 and 
8 representing typical values of 5 q and 8 q . The third term in Eq. (I]) [25] is the 
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gCEDM, with 



in terms of a dimensionless parameter w. The sizes of 5, 5 and u> depend on the 
exact mechanisms of electroweak and T breaking and on the running to the low 
energies where non-perturbative QCD effects take over. The minimal assumption 
is that they are 0(1), 0(g a /4:ir) and (D((g s /4:ii) 3 ), respectively, with g s the strong- 
coupling constant. However they can be much smaller (when parameters encoding 
T-violating beyond the Standard Model are small) or much larger (since f q is unnat- 
urally small). In the Standard Model itself, where Mr = M w , they are suppressed [H] 
by the Jarlskog parameter [29J Jqp — 3 • 10~ 5 . In supersymmetric models with var- 
ious simplifying, universality assumptions of a soft-breaking sector with a common 
scale Msusy, one has Mf = Msusy and the size of the dimensionless parameters 
is given by the minimal assumption times a factor which is [21 [30l El] , roughly (ne- 
glecting electroweak parameters), Aqp = (g s /4:7i) 2 sin</>, with <p a phase encoding 
T violation. Allowing for non-diagonal terms in the soft-breaking sfermion mass 
matrices, enhancements of the type m^/md ~ 10 3 or even m t /m u ~ 10 5 are possible 
(although they might be associated with other, smaller phases) [2]. 

Since we are interested in light systems, we are integrating out all degrees of 
freedom associated with quarks heavier than up and down. The effects of qEDMs 
and qCEDMs of such quarks are discussed briefly at the end. T-violating four-quark 
operators are effectively dimension-8 because again electroweak gauge invariance re- 
quires insertions of the Higgs field. Since higher-dimension operators are suppressed 
by more inverse powers of the large scale Mr, we expect them to be generically less 
important at low energies and we concentrate here on the dimension-6 operators in 
Eq. Q. It is of course possible that in particular models the coefficients of the 
effectively dimension-6 operators are suppressed enough to make higher-dimension 
operators numerically important; low-energy implications of four-quark operators, 
which also contain representations of chiral symmetry we consider, have recently 
been studied in Ref. [32] . 

At momenta Q comparable to the pion mass, Q ~ m n <C Mqcd, interactions 
among nucleons, pions and photons are described by the most general Lagrangian 
involving these degrees that transforms properly under the symmetries of the QCD. 
Ignoring quark masses and charges and the 9 term, the dimension-4 QCD terms 
are invariant under a chiral SU{2)l x SU(2)h ~ 5*0(4) symmetry, which is spon- 
taneously broken down to its diagonal, isospin subgroup SU(2) V ~ SO(3). The 
corresponding Goldstone bosons are identified as the pions, which provide a non- 
linear realization of chiral symmetry. Pion interactions proceed through a covariant 
derivative, which in stereographic coordinates p2] w for the pions is written as 

D^n = ZT^tt, (4) 

with D = 1 + ir 2 /F% and F n ~ 186 MeV the pion decay constant. Nucleons are 
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described by an isospin-1/2 field N, and the nucleon covariant derivative is 



We define X>t through NV^ = VN, and use the shorthand notation 
and 

Covariant derivatives of covariant derivatives can be constructed similarly, for ex- 
ample 

(ZV^tt); = (^<% " ^^'(w x D ^)k) D ^y (8) 

For simplicity we omit the delta isobar here, although one can introduce [33] an 
isospin-3/2 field for it along completely analogous lines. The effective interactions 
are constructed as isospin-invariant combinations of chiral-covariant objects [T7] . 

The quark mass, charge and 9 terms break chiral symmetry explicitly as specific 
components of various chiral tensors. The formalism to include chiral-symmetry- 
breaking operators in the SU(2) x SU (2) ChPT Lagrangian has been developed in 
Refs. [171 [33]. Introducing the SO (A) vectors 

s= ( -%qi°Tq \ p= ( qrq \ , Q 

\ qq J ' \ iqrfq J ' 

and the SO (4) scalar and antisymmetric tensor 

1 - j* n ™ _ 1 / ^jkqYl 5 r k q q^Tjq 



p = -q^ q , T M = - ^ fc<i _' ' y ' , (10) 

6 2 V -qi^m J 

the average quark-mass term transforms as S4, the quark-mass-difference term as 
P 3 , the quark-photon coupling as / © T 34 , and the 9 term as P 4 . They generate 
interactions containing the pion field explicitly, which are proportional to powers 
of the symmetry-breaking parameters fh = (m u + m<i)/2, em = {ma — m u )/2, e 
(the proton charge), and (m(l — £ 2 )sin#)/2. The most important chiral-breaking 
term is the m term, which among other effective interactions generates the main 
contribution to the pion mass, m\ = (9(?7iMqcd)- The electromagnetic coupling 
produces two types of effective interactions: i) purely hadronic interactions propor- 
tional to a em /47r ~ em^/MQ CD ; and ii) gauge-invariant interactions with explicit 
soft photon fields, which appear either in gauge-covariant derivatives or through the 
photon field strength. The covariant derivatives below are all to be interpreted as 
gauge-covariant derivatives. After a suitable chiral rotation eliminates it in favor of a 
mass term that does not generate vacuum instability in first order in the symmetry- 
breaking parameter [31] , the 9 term is found to break chiral symmetry as a different 
component of the same vector P to which the isospin-breaking quark mass term is 
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associated. The construction of the corresponding effective interactions has been 
carried out in some detail recently [TH]. Since effective interactions proportional 
to two or more powers of T-violating parameters are exceedingly small, to a very 
good approximation one can simply add the contributions from dimension-6 sources 
considered here to the corresponding 6 contributions calculated in Refs. [TO} [12]. 

Since nucleons are essentially nonrelativistic for Q <C m^, the nucleon mass, we 
work in the heavy-baryon framework [16J where, instead of gamma matrices, it is 
the nucleon velocity v M and spin 5 /t (5 = (<?/2, 0) in the rest frame v = (0, 1)) that 
appear in interactions. Below we use a subscript _L to denote the component of a 
four-vector perpendicular to the velocity, for example 

V\ = XX 1 - v^v ■ V. (11) 

We use reparametrization invariance (RPI) [35] to incorporate Lorentz invariance in 
an expansion in powers of Q/m N . 

The infinite number of effective Lagrangian terms can be grouped into sets 
of a given "chiral index" [15] A = d + f /2 — 2, where d counts derivatives, powers 
of m n and photon fields, and / = 0, 2 is the number of fermion fields: 

oo 

£ = ^£( A ). (12) 
A=0 

The LECs can be estimated using naive dimensional analysis (NDA) [251 ES], in 
which case the index A tracks the number of inverse powers of Mqcd ~ 27ri^ ~ 1.2 
GeV associated with an interaction. (Note that since NDA associates the LECs of 
chiral-invariant operators to g s /4:7c, for consistency one should take g s ~ Air.) For 
the purposes of our calculation, we need explicitly only the leading T-conserving 
interactions, 

£ {0) = ■ D»n - ^tt 2 +N(ivV- ■ D^j N, (13) 

where — 1.267 is the pion- nucleon axial- vector coupling. At this order the nu- 
cleon is static; kinetic corrections have relative size 0(Q/Mqqd) and appear in CP-' . 
Isospin breaking from the quark masses, represented by e, also only appears in 
subleading orders [33] . 

The dimension-6 sources of T violation generate further effective interactions, 
which break chiral symmetry in their own ways. Introducing the 50(4) singlet 

Iw = \^ Xo ' f ahc Gl p G b /G c Xa , (14) 

and the 50(4) vectors 

W -2{ qa^q ) F ^ *\ iqa^q ) F ^ (15) 

and 
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Eq. ([I]) can be rewritten as 

C f = -d 1/ 4 + d 3 W 3 - d V 4 + d 3 W 3 + d w I w . (17) 

The corresponding T- violating effective Lagrangian can be constructed by writing 
down all terms that transform in the same way under Lorentz, P, T, and chiral 



symmetry as the terms in Eq. (17). We still use NDA and label operators by a 
generalized chiral index A that continues to count inverse powers of Mqcd- For 
simplicity we will not keep track of explicit factors of e. Here we present only the 
interactions needed in the calculation of the nucleon EDFF up to the order the 
SM first appears; as we are going to see, this means up to A = 1 for qCEDM 
and gCEDM and to A = 3 for qEDM. In the equations below, ". . . " account for 
interactions not needed in our calculation; we leave a more complete presentation 
of the effective Lagrangian for a future publication [28] . 

Some of the contributions to the EDFF arise from virtual pions. In the presence 
of T violation, pions can be annihilated into the vacuum when an operator with 
the quantum numbers of the neutral pion is allowed by the pattern of symmetry 
breaking. In the case of the qCEDM, such a tadpole arises from W 3 and can thus be 
linear in d 3 . In the case of the gCEDM, the tadpole arises from the tensor product of 
Iw with the P 3 in Eq. (|9| and is linear in efndw- In both cases these tadpoles first 
appear at A = — 2 and exist also at higher orders. For the qEDM, they are beyond 
the order we consider here because they are suppressed by at least one factor of 
a e m/47r. All such tadpoles represent vacuum misalignment. Because they are small, 
they can be treated in perturbation theory or simply eliminated using the chiral 
rotation given in Ref. [19J. To the order we are working the effects of this rotation 
can be absorbed in terms that already exist in the effective Lagrangian. The fields 
and LECs introduced below are to be interpreted as subsequent to the rotation. 

Pions contribute to the EDFF in loops, which require the T-violating pion- 
nucleon interactions with A = — 1. Again, "indirect" electromagnetic operators 
stemming from hard photons tied to qEDM are of higher order. From the qCEDM 
and the gCEDM, 

47^ = -J^k ■ NtN -^(v-Dttx D^) ■ NS^tN + ... (18) 



The non-derivative term in Eq. ( 18 ) arises either directly from V4 or from the tensor 



product Iw ® 5*4, and thus has a LEC 

^M QCD! u^M QCD | . (19) 

The two- derivative term is the lowest-index chiral invariants that arise from Iyy, its 
LEC being 

*o = ^MqcdJ . (20) 

Note that the qCEDM and the qEDM generate two- derivative interactions of differ- 
ent form than above, since they are chiral-breaking, but they only appear at higher 
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order. There are other pion-nucleon interactions with A = — 1, but they do not 
contribute to the EDFF at the order we calculate. 

In addition to the long-range contributions from virtual pions, the EDFF is sensi- 
tive to shorter-range effects, which in ChPT are represented by contact interactions. 
The lowest-order contribution of this type arises from the gCEDM combined with 
the quark-photon coupling, Iw © (I © T34): 



where the LECs are 



T3 + ^7j (*-37r-T-7T 2 T 3 ) 



S"Nv v F^ + ..., (21) 



= O 



ew 



Mr 



2^QCD 



T 



In next order, there is a recoil correction 



C, 



(0) _ 1 ^-^^n(- 1 ) , n(-i). 



m N 



N D K ' + D [ f >r 3 S»V V ± _N Fp V + . . . , 



and one further order up other sources contribute as well: 



-(i) 



= 2N 



D, 



2tt 2 \ 



(i) 



2tt, 



T 3 



FID 



77 • T 



S^Nv"F, 



-2D\ 



(i)' 



2tt 2 \ - 



T 3 



FID 



(7T 3 7r • T - 7T 2 T 3 ) 



(22) 



(23) 



r w 



-N (s} x) + ^ (1 V 3 ) {S ■ V ±+ V» ±+ + S»V 2 ±+ ) Nv»F, u 



4m% 



N ( L>S _1) + D\ l 't 3 ) S ■ V^V^_Nv v F^ v + ... 



(24) 



Here the operators with LECs 



Mr 



T 



Mr 



T 



(25) 



transform as vectors: the isoscalar component as V4 or as the vectors in © I, 
W3 © T34 and Iw © (£4 © P 3 ) © (/ © T 34 ); the isovector component as W3 or as the 
vectors in W 3 © I, V 4 © T 34 and V © (£4 © P 3 ) © (/ © T 34 ). The operator with LEC 



(i)' 



(26) 



transforms as the tensors in V 4 © T 34 and 7w © S 4 © (7 © ^34) ■ The operators with 
LECs 



S'P = O 



ew x 

M} Q° d 



(27) 
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come from Iy/ <g> (J © T 34 ). The last operator written explicitly in Eq. (24) is a 



relativistic correction. It is important to realize that the form of such corrections 
depends on the choice of operators included in the effective Lagrangian; here we 
have eliminated time derivatives of the nucleon field through field redefinitions. 
For the qEDM, we need also 



A*) 



m N 



-N 



ST> v ] _NF tH , + ... 



and 



A3) 



with 



2N (B® + D^n} S^NvT^ 

-N (£? 3) + S; (3) r 3 ) (S ■ V ±+ V<* ±+ + S»V\ + ) Nv u F f 



(2f 



Am 2 N 



N ( + D^t 3 ) S ■ V JL JD*__Nv v F ia , + 



(29) 



D 



(3) 







TT) 



5^ 



/(3) 



(30) 



With these interactions we can calculate the nucleon EDFF to the order at which 
momentum dependence first appears. We consider a nucleon of initial (final) mo- 
mentum p (p') and a (space-like) photon of momentum q = p — p' (q 2 = —Q 2 < 0). 
It is convenient to take q and K = (p + p')/2 as the independent momenta. The 
isoscalar (Fq) and isovector (F\) EDFFs are defined from the nucleon electromag- 
netic current Jg m (q) via 

Jgjq, K) = 2 (F (Q 2 ) + i^QVs) {S"v ■ q - S ■ qv» 

1 1 1 

+ [q ■ KS» - S ■ qK") + ^-yS • K [K»v -q-K- qv") +... \(31) 



m N 



2m 



N 



The first term corresponds to the definition in Ref . [TU] , while the second is a recoil 
correction [TS] and the remaining are consequences of Lorentz invariance. We will 
write 



F l (Q 2 ) = D t -S' i Q 2 + H t (Q 2 ), 



(32) 



where Di is the isospin i component of the EDM, £• of the SM, and Hi(Q 2 ) accounts 
for the remaining Q 2 dependence. The EDFF of the proton (neutron) is Fq + F\ 
(F -F 1 ). 

The calculation of the EDFF to the order we are interested in includes T violation 
in tree and one-loop diagrams. In tree diagrams the photon is attached to the nucleon 



line via a T- violating operator from Eqs. (21), (23), (24), (28), and (29). The loop 



diagrams, shown in Fig. [TJ contain the T-violating pion-nucleon couplings in Eq. 
(18) or the photon-nucleon couplings in Eqs. (21) and (24) — which we denote 



by squares — while the other couplings come from the leading, T-preserving chiral 



Lagrangian, Eq. (13). In addition, nucleon wave-function renormalization from 
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<vw> ■ WV 



5 



7 



IK/W> 

8 
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FIG. 1: One-loop diagrams contributing to the nucleon EDFF. Solid, dashed and wavy 
lines represent the propagation of nucleons, pions and photons, respectively. A square 
marks a T-violating interaction, other vertices representing T-conserving interactions. 



at one-loop level [IE] can appear. Since in this Lagrangian the nucleon is static, in 
one-loop diagrams we take v ■ q = v ■ K = 0. We use dimensional regularization in 
d dimensions and encode divergences in the factor 



4-d 



— 7s + ln47r . 



(33) 



The loops bring in also a renormalization scale //, which is eliminated through the 
accompanying LECs. 

We start with the contributions from the qCEDM, which are very similar to those 
of 9 [TUl H2] . In this case the lowest-order momentum dependence arises from the 
loops 1-6 in Fig. [TJ where the T-violating vertex is the first term in Eq. (18). At 
the same order, 0(e5m 2 /MpMQCT>), there are also tree contributions from the first 
term in Eq. (24). The isoscalar form factor does not receive loop corrections and 



can be expressed purely in terms of coefficients of short-distance operators, 



A),qCEDM = A 
^0,qCEDM = 0; 
-^0,qCEDM(Q 2 ) = 0. 



(1) 



(34) 
(35) 
(36) 



Contributions to the isoscalar SM appear in higher orders. In contrast, the loop dia- 
grams with static nucleons not only renormalize the contributions of short-distance 
operators to the isovector EDM, but also generate a non-trivial momentum depen- 
dence in the isovector EDFF. The /i-independent isovector EDM is found to be 

,2 



-t^l.qCEDM — L>\ + U \ + 



(27TF ff )2 

while the momentum dependence is encoded in 



L-ln 



q/ 

°l,qCEDM 
-f^l,qCEDM(Q 2 ) 



6774(27^)2 



15(2^)2 



./ 



where the function f(Q 2 /4m^) is defined as 



15 

T 



In 



x 



4777,2 / ' 



-2 1 + 



(37) 

(38) 
(39) 



(40) 
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Note that f(x < 1) = x 2 + 0(x 3 ). 

Contrary to the qCEDM, the momentum dependence for qEDM and gCEDM 
arises only two orders down with respect to the lowest-order contribution to the 
EDM. To this order, a calculation of the electromagnetic current yields, in addition 
to strong-interaction corrections, also the Lorentz terms oc and oc m~^ 2 in Eq. 



(31 ). In the strong-interaction corrections given below we include the nucleon wave- 



function renormalization. 

For the qEDM short-range contributions to the EDM start at chiral index A = 
1 and others appear at A = 3. To this order there are no contributions from 
pion-nucleon T-violating interactions, while the loop diagrams 7, 8, 10, and 11 in 
Fig. [TJ with interactions from Eq. (24), only renormalize the tree-level contributions 
without any energy dependence. To 0(e<Jm*/M|Af$ CD ), we find the EDMs 



-Do, qEDM 
-Dl >q EDM 



r)(l) I r)W I - 



,(3) 



mz 



j (3) 



mz 



(2 + 4^) L-ln-f )+2 + g A 



mz 



(2 + 8^1) (L-ln-f ) 2 + 3,/-, 



(41) 
,(42) 



and the momentum dependence given entirely by the SMs, 



q> 

°j,qEDM ' 
//. 



s 



/(3) 



'i.qEDM 



(Q 



o. 



(43) 
(44) 



In the case of the gCEDM, short-range contributions to the EDFF start at A = 
— 1, which dominate, and appear again at A = 1, suppressed by m 2 /MQ CD . At this 
order there are also contributions from the T-violating pion-nucleon interactions in 
Eq. (18) through the loops 1-10 in Fig. [TJ and from the photon-nucleon interactions 



in Eq. (21) through the loops 10 and 11. Thus, to 0(ewm 2 /MZMQCT>) we find the 
/^-independent EDMs 



-Do,gCEDM 
-Dl.gCEDM 



D 



(-1) 



mz 



(2ttF 7T 



mz 



+ 



{l+2g 2 A ) D[- 1] + e 



V mi 

\ 7T 



L-ln 



mz 



(2<kF v ) 
( Ma , io 



+ 



ei 



L-ln 



mz 



(45) 



(46) 



The isoscalar momentum dependence is entirely due to short-range operators in Eq. 



(24), 



ql 

°0,gCEDM 



-^O.gCEDM (Q ) = o. 



(47) 
(48) 
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The isovector part, on the other hand, receives also non-analytic contributions: 



We are now in position to discuss the implications of the various dimension-6 
T- violation sources to the nucleon EDFF. 

First, we note that to this order the nucleon EDFF stemming from the qCEDM 
has a form that is identical to that [TU1 [12] from the 9 term. In both cases the 
momentum dependence (and thus the SM) is isovector, has a scale (relative to the 
EDM) set by Im^, and is determined by the lowest-order pion nucleon coupling g . 
The EDFF depends on just three independent combinations of LECs, go and the 
short-range EDM contributions and + , which contain nucleon matrix 
elements of V4 for qCEDM and P4 for the 9 term. The numerical factors relating 
these couplings to either 5 or 9 will thus be different. In the case of 9, the matrix 
element in g can be determined from T-conserving observables, because it is related 
[T9"] to the matrix element of P3 that generates the quark-mass contribution to the 
nucleon mass splitting: g /9 ~ 3 MeV. For the qCEDM, an argument identical 
to that in Ref. [20] serves to estimate -Di^cedm in terms of g , but no analogous 
constraint exists for go in this case and without a lattice calculation or a model we 
cannot do better than dimensional analysis. (For an estimate with QCD sum rules, 
see Ref. [39].) In any case, to the order we consider here, any EDFF measurement 
alone will be equally well reproduced by a certain value of 9 or a certain value of 5. 
Note that the qCEDM does give rise to additional effective interactions generated by 
W3, which contribute to the nucleon EDFF only at higher orders but could generate 
sizable differences for other observables. 

Second, the pion-nucleon sector of the qEDM is suppressed compared to that of 
the qCEDM because of the smallness of a em compared to g 2 /^ at low energies. The 
consequence is that, up to the lowest order where momentum dependence appears, 
both the EDM and the SM from the qEDM are determined by four combinations 
of six independent LECs, which at this point can only be estimated by dimensional 
analysis. The momentum dependence is expected to be governed by the QCD scale 
Mqcd, small relative to the EDM, and nearly linear in Q 2 . 

Finally, in the case of the gCEDM loops are also suppressed, but do bring in 
non-analytic terms not only to isoscalar and isovector EDMs, but also to the isovec- 
tor momentum dependence (and thus SM). Again the momentum dependence is 
governed by Mqcd- in addition to seven short-range contributions to the EDMs 
and SMs, also two independent pion-nucleon LECs appear (g and z ) which endow 
the isovector EDFF with a richer momentum dependence than in other cases. The 
isoscalar momentum dependence is identical to qEDM. For the gCEDM, using the 
pion loop together with an estimate of go [40J is likely to be an underestimate of the 




(49) 



(50) 
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TABLE I: Expected orders of magnitude for the neutron EDM (in units of c/Mqcd); the 
ratio of proton-to-neutron EDMs, and the ratios of the proton and isoscalar SMs (in units 
of 1/m 2 ) to the neutron EDM, for the 9 term [10, 12J and for the three dimension-6 sources 
of T violation discussed in the text. 



Source 



MqcDd n /e 
dp/d n 
mlS' p /d n 

mlS' /d n 



qCEDM qEDM 



gCEDM 



O 9 




O 



w 



M QCD 



Oil) 



o 



o 



M. 



QCD 



EDM, because chiral symmetry allows a short-range contribution that is larger by 
a factor M 2 CD /ml 



As it is clear from Eqs. (|36j), (|39j), (|44j), (|48j), and ([SfJ, the full EDFF momentum 
dependences (for example, the second derivatives of Fi with respect to Q 2 ) are differ- 
ent for qCEDM (and 6), qEDM, and gCEDM. Although the isoscalar components 
all have linear dependences in Q 2 (with different slopes) to the order considered 
here, the isovector components show an increasingly complex structure as one goes 
from qEDM to 6 and qCEDM to gCEDM. Determination of nucleon EDMs and SMs 
alone would not be enough to separate the four sources, yet they would yield clues. 
Expectations about the orders of magnitude of various dimensionless quantities are 
summarized in Table IB 

In the first line of Table|T]one finds the expected NDA size of the neutron EDM. As 
it is well known [2], this is consistent with many other estimates, such as d n = 0{dj) 
in the constituent quark model, and d n = O^edi/iii, edw A^qcd / '47r) from QCD sum 
rules. If 5 ~ 5 ~ w = 0(1) (as would be the case for g s ~ Att and no small 
phases), then the gCEDM gives the biggest dimension-6 contribution to the EDFF 
because of the chiral-symmetry-breaking suppression 0(m 2 /MQ CD ) for the qCEDM 

and qEDM. However, models exist (for example, Ref. [30]) where 5 and 5 are 
enhanced relative to w, and all three sources produce EDFF contributions of the 
same overall magnitude. Even so, there is no a priori reason to expect cancellations 
among the various sources. A measurement of the neutron EDM d n could be fitted 
by any one source. Conversely, barring unlikely cancellations, the current bound 
yields order-of-magnitude bounds on the various parameters at the scale where NDA 
applies: using 2TrF n ~ 1.2 GeV for Mqcd, 

< 10- 10 , (51) 



5 5 
M^'M 2 
w 
M 2 



< (l0 5 GeV) 2 , (52) 

< (l0 6 GeV)~ 2 . (53) 



(For comparison, Eq. (51 ) is consistent within a factor of a few with bounds obtained 
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by taking representative values of \x in the non-analytic terms to estimate [20] the 
size of the renormalized LECs for the EDM, and using either SU (2) [TH] or 577(3) 
[12] symmetry to constrain go.) In all four cases we expect the proton and neutron 
EDMs to be comparable, \d p \ ~ \d n \, but the presence of undetermined LECs does 
not allow further model-independent statements. 

It is in the pattern of the £• that we see some texture. (This pattern is not evident 
in Ref. [26], possibly because of the way chiral symmetry is broken explicitly in the 
model used, both in the form of the T-conserving pion-nucleon Lagrangian and 
in the gCEDM magnitude of the T- violating pion-nucleon coupling.) While in all 
cases one expects \S' p \ ~ |S^|, the relative size to the EDMs, in particular of the 
isovector component, allows one in principle to separate qEDM and gCEDM from 
9 and qCEDM. Since all these sources generate different pion-nucleon interactions 
thanks to their different chiral-symmetry-breaking properties, nuclear EDMs might 
provide further probes of the hadronic source of T violation. 

More could be said with input from lattice QCD. For each source the pion-mass 
dependence is different. A fit to lattice data on the Q 2 and ml dependences of 
the nucleon EDFF with the expressions of this paper would allow in principle the 
separate determination of LECs. In this case a measurement of the neutron and 
proton alone would suffice to pinpoint a dominant source if it exists, but in the more 
general case of two or more comparable sources further observables are needed. 

One should keep in mind that our approach is limited to low energies. The 
contributions associated with quarks heavier than up and down are buried in the 
LECs, as done, for example, in other calculations of nucleon form factors: electric 
and magnetic [37], anapole [38], and electric dipole from 9 [TO] . Heavy-quark EDMs 
and CEDMs are also singlets under SU(2) L x SU(2) R , so they generate in two-flavor 
ChPT interactions with the same structure as those from the gCEDM, and cannot 
be separated explicitly from the latter. (This is clear already in the one-loop running 
of d\y, which gets a contribution of the heavy-quark CEDMs [22] ■) The parameter 
w here should be interpreted as subsuming heavier-quark EDMs and CEDMs. With 
the additional assumption that m s makes a good expansion parameter, effects of 
the s quark could be included explicitly. The larger S77(3)x x SU(3)r symmetry 
would yield further relations among observables (for example, between the EDFFs 
of the nucleon and of the A), and we could, in principle, isolate the contributions of 
the strange quark. Since our nucleon results, which can be used as input in nuclear 
calculations in two-flavor nuclear EFT, would be recovered in the low-energy limit 
anyway — as was explicitly verified in Ref. [T2] for the 9 results of Ref. [JU] — we 
leave a study of the identification of explicit s-quark effects to future work. 

In summary, we have investigated the low-energy electric dipole form factor that 
emerges as a consequence of effectively dimension-6 sources of T violation at the 
quark-gluon level: the quark electric and color-electric dipole moments, and the 
gluon color-electric dipole moment. Only the full momentum dependence could in 
principle separate these sources, although the Schiff moments, if they were isolated, 
would partially exhibit a texture of T violation. Further implications of the different 
chiral-symmetry-breaking patterns of these sources will be studied in a forthcoming 
paper [281 . 
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